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Abstract. Lower bounds estimates are proved for the first eigenvalue for the 
Dirichlet Laplacian on arbitrary triangles using various symmetrization tech- 
niques. These results can viewed as a generalization of Polya's isoperimetric 
bounds. It is also shown that amongst triangles, the equilateral triangle min- 
imizes the spectral gap and (under additional assumption) the ratio of the 
first two eigenvalues. This last result resembles the Payne-Polya- Weinberger 
conjecture proved by Ashbaugh and Benguria. 



1. Introduction 

The purpose of this paper is to prove new isoperimetric-type bounds for the eigen- 
values of the Dirichlet Laplacian on arbitrary triangles. Given a domain D we 
will use "eigenvalue of the domain D" to refer to an eigenvalue of the Dirichlet 
Laplacian on the domain D. The eigenvalues for a domain form a nondecreasing 
sequence {Xi,j > 1} with Ai < A2. Throughout the paper we use A for the area 
of a domain, L for its perimeter, R for its inradius (the supremum of the radii of 
disks contained in the domain) and d for its diameter. 

The problem of finding good bounds for eigenvalues of various domains has been 
of interest for many years. See for example [1,^9, ^15^^16^ for general bounds, |2 „2T] 



for result about polygons, and 10 11 18 for bounds for triangles. Comprehensive 



overview of this subject along with methods used to tackle it can be found in the 
book The results contained in [11 are especially interesting, since they are 

of different nature than our bounds and it should be possible to combine methods 
from this paper with our approach to get even better bounds. 



In 15 Section 7.4], Polya and Szego conjectured that with fixed area, amongst 
all polygons with n sides the regular one minimizes the first eigenvalue (they also 
conjectured that the inner radius and the transfinite diameter are also minimized, 
but these has been proved). The conjecture remains open except for the following 
cases: 

(1.1) XlA\T„angle > XlA\ 

Equilateral i 

(1.2) \iA\ 

Quadrilateral Square •> 

(1.3) \lA>\^A\BaU- 

The notation \set is used to indicate the set for the quantity to the left of it. In the 
last bound (called the Faber-Krahn inequality) no set is specified since the bound is 
true for an arbitrary domain. The proofs of these results along with the conjecture 
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for polygons can be found in ^12j. Note that the ball in the last bound can be 
viewed as a limiting case of a regular polygon with an infinite number of sides. 
There are also upper bounds where in this case the ball is the extremal. 



(1.4) 
(1.5) 

(1.6) 



Ball, 



Ai 



< 



Ai 



Ball 



{\2-\i)B? < (A2-Ai)i?2| 



Ball- 



The first inequality follows trivially from the domain monotonicity. The second and 
the third are known as the Payne-Polya- Weinberger conjecture proved by Ashbaugh 
and Benguria T]. Note that there is no scaling factor (area or inradius) in the second 
bound. The difference of the eigenvalues in the last bound is called the spectral 
gap and it is important in the study of dynamical systems. It can be regarded as a 
measure of the speed of convergence to equilibrium. For more on this, see [s] and 
19]. 

We can see that a ball gives both upper and lower bounds for general domains. 
We want to show that an equilateral triangle has exactly the same properties among 
triangles. We already have the lower bound ( |1.1[ ) as an analog of (1.3 1. It is proved 
in (Tsl that 



(1.7) 



Aii?2 



\Triangle 



uilateral • 



This is the parallel of the bound (1.4 1. The first result of this paper is the following 
Theorem 1.1. 

Ai)i?2| Equilateral • 



For arbitrary triangle 

(1.8) (A2 — Xi)R^\Triangle < (A2 

// a triangle is acute we also have 



(1.9) 



A2 
Ai 



< 



A2 
Ai 



Equilateral 



Triangle 

The additional assumption in the second part of the theorem is due to the method 
used to prove the result. We believe this result should hold for all triangles. The 
proof of the theorem relies on a variational formula for eigenvalues, as well as some 
very cumbersome computations. 

In view of these results, we venture to propose the following generalization of 
Polya's conjecture about polygons 



Conjecture 1.2. 

gon with n sides. 

(1.10) 
(1.11) 
(1.12) 

(1.13) 



Let P{n) denote a polygon with n sides and R{n) a regular poly- 
Then 



(A2 



Ai^|p(„) > AiA|j^(„), 

XlR?\p{n) < Aii?^ 
Al)^^|p(n) < (A2 - 



Ai)i?2| 



i?(n)j 



A2 

Ai 



< 



P{r 



A2 

Ai 



R(r 



It is easy to check that this conjecture is true for rectangles and that for quadri- 
laterals (1.10) is the same as (1.2). Our results along with ( |1.1[ ) and (1.7) prove 
this conjecture for triangles, except for obtuse triangles in the last bound. All other 
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cases remain open. It is worth noting that a shghtly weaker version of is a 

part of 21 Theorem 2]. The only difference is that the scahng factor R is missing 



and the regular polygon has outer radius 1 instead of inradius 1 . As we have been 
recently informed by the author of |21|, the proof of this theorem should work if we 



replace outer radius 1 with inradius 1. This would prove inequality (1.11 1 



The second goal of this paper is to establish sharper lower bounds for the first 
eigenvalues of arbitrary triangles. In addition to ( |1.1[ ) we have 

^2 



(1.14) 
(1.15) 



Al > 7T^ 



4 



4/12 / 



The first bound is due to Protter 
Freitas 
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16 and the second result was proved recently by 
Each of these bounds is the better than the other bounds for certain 



triangles, but not as good for some others. 



We obtained new lower bounds for triangles which are better than (1.141 and 



(1.151 whenever these are better than (1.1) 



Theorem 1.3. For an arbitrary triangle with area A, diameter d and shortest 
altitude h, 



(1.16) 



d2 + 



4A2 



We also obtained a sharp bound based on circular sectors. 

Theorem 1.4. Let 7 be the smallest angle of a triangle. Denote by /(A, 7) an 

isosceles triangle with same area A and the vertex angle 7, and S{A, 7) a circular 
sector with area A and angle 7. Then 

(1-17) \l\Triangle > Ai|/(^^^) > Xl\s{A,j)- 

The function 

(1-18) f{l) = M\liA,y) 

is decreasing for 7 G (0, 7r/3) and increasing for 7 G (tt/S, tt). 

This last result can be viewed as a generalization of the Polya's isoperimetric 
bound (1.1). If we fix A and the smallest angle, then the isosceles triangle min- 
imizes the first eigenvalue. Then, due to the monotonicity property, we also get 
an alternative prove of ( |1.1[ ). The bounds involving circular sectors can be used to 
get good lower bounds for the eigenvalues since the eigenvalues of sectors are given 
explicitly in terms of the zeros of the Bessel function. 

From now on I{A, 7) will always denote an isosceles triangle with area A and 
angle 7 between its sides of equal length. This angle will be called the vertex angle. 
The side opposite to that angle will be called the base and the other two the arms. 



The methods used to prove (1.15) and (1.14) are not based on any kind of 
symmetrization argument. In contrast, the proofs of our lower bounds rely on 
certain symmetrization techniques similar to those used in the proof of Polya's 
isoperimetric inequality. In this sense, our results can be viewed as generalized 
isoperimetric bounds. Using the same techniques we also give an alternative proof 
of Freitas's bound (1.15). This shows that symmetrization is the best way for 
obtaining lower bounds for the eigenvalues, at least for triangles. The bound is 
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determined by the shape of the symmetrized domains. That is, the equilateral 
triangle in (l.ll, the rectangles in Theorem 1.3 and (1.151, or the circular sector in 



Theorem 11.41 

The rest of the paper is organized as follows. In Section 2 we compare the new 
lower bounds with the known results. Section 3 contains the definitions and ba- 
sic explanations of various forms of symmetrization techniques which are used in 
Section 4 to prove the lower bounds. Section 5 contains the proofs of the upper 
bounds. Variational methods are used to obtain complicated polynomial bounds 
for the eigenvalues. These bounds are then simplified to the final results by proving 
certain polynomial inequalities. The algorithm for solving such polynomial inequal- 
ities is given in Section 6. The last section contains a script in Mathematica used 
to perform long calculations. 



2. Comparison of lower bounds 



In this section we show that our bound (Theorem 1.3 1 is sharper than the bound 
obtained by Freitas (1.151 whenever the latter is better than Polya's isoperimetric 
bound (1.1). We also give a numerical comparison between the lower bounds to 
show that Theorem 1.4 in practice gives the best lower bounds for a wide class of 
triangles. 



Both bounds (1.15) and Theorem 1.3 are of the form 

'4 X \ 



(2.1) fix) = TT^ 

If we write the explicit value for the eigenvalue of the equilateral triangle, the bound 
(l.ll reads 



(2.2) 



Ai > 



4^3^ 
3A 



To compare it to the other two we need to investigate the following inequality 
(2.3) 

Put X = AyA, then 
(2.4) 



1^ > 1 
3A - X 



X 

4A2' 



4^3 ^ 1 
3 - y 



y 
1' 



One ca n ch eck that the equality holds if j/ = VS o r y = \/3/3. Hence, Polya's 
bound (O) is better than a bound of the type (O) if x G {AAV3/3,AAV3). We 
also observe that f{x) is increasi ng for x > AAy/S. 

In the case of Freitas's bound (1.15) (proved in jio]) we get x — d? . Hence this 
bound is better than Polya's bound if 



(2.5) 



^ (^4ylV3/3, 4AV3) 



If we denote the length of the altitude perpendicular to the side of length d by /i 
we see that 



(2.6) 



(2/i\/3/3,2/i\/3 
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Observe that d cannot be smaller than 2ft,\/3/3 and it is equal to this quantity only 
for an equilateral triangle. Hence Freitas's bound is better than Polya's bound if 

(2.7) 



In Theorem 1.3 we have x 
bound 



d > 2hV3- 

f + which is a bigger value than in Freitas's 
Therefore this bound is the best for every triangle such that the above 
condition is true. Figure [T] sh ows where each of the lower bounds (Protter ( |1.14 1, 
Polya (1.1 1 and Theorem 1.3 1 is the best. 



Theorem 1.3 





M 3 



Figure 1: Theorem 1.3 



On this, and all other figures, M denotes the side with the middle length, and 
U is equal to the difference between the length of the longest side and M. The 
shortest side is assumed to be 1, hence U < 1. This gives a one-to-one mapping of 
all the triangles onto the infinite strip [0, 1) x [l,oo). 

Next, we give some numerical results with Theorem |1.4| included. Since the 
bound involving the eigenvalues of circular sectors rely on the calculations of the 
zeros of the Bessel function, it is hard to compare to the other bounds. The nu- 
merical comparisons from Figure 2a show that this bound is better than the other 
bounds for almost all triangles. 

One could also use a simplified sector based lower bound by just taking the 
smallest sector containing a given triangle. This gives the sector S{'yN'^ /2,j) with 
the same angle as in Theorem 1.4 but larger radius. This is clearly worse than 
Theorem lA but still gives a good bound as can be seen on Figure |2b] 

We do note here that there is no best bound, although our new bounds together 
with Polya's give the best results depending on the triangle. We also see that 
symmetrization techniques (all bounds are consequences of certain symmetrization 
procedures) are very powerful and lead to very good lower bounds for the first 
eigenvalue of triangles. 

It would be also interesting to see how far the bounds are from the exact values. 
We begin with a right isosceles triangle with arms of length 1. The first eigenvalue 
is known in this case and it is equal to Stt^. The lower bounds results are given in 
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1 Af 3 5 

(b) Smallest sector containing a triangle 



Figure 2: Sector-based bounds 



Table [Ta] Clearly, neither bound is particularly accurate, although Polya's bound 
and Theorem 11.41 are the closest. 

The latter bound works best for "taller" triangles. Consider a right triangle with 
angle 7r/6 and hypotenuse 2. Table lb shows the values of the lower bounds for 
this triangle. This time Theorem |1.4| gives a very close bound, while the others 
are not so accurate. The advantage of the sector bound is the biggest for very 
"tall" triangles. For an isosceles triangle with base 1 and fixed arms we get Table 
[2] The exact values are not known, but the difference between the bounds given 
by Theorem 1.4 and the other bounds is clear. 

All these numerical results show that the lower bounds can be very accurate if 
the triangle is acute and almost isosceles (Polya's bound or Theorem 1.4). Unfor- 
tunately, neither bound is very good in the case of "wide" obtuse triangles. The 
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exact 

Polya (|1 J| 
Freitas (05 1 
Protter ( 
Theorem 
Theorem 



49.348 

45.5858 

39.4784 

29.9958 

40.4654 

45.2255 



(a) Right isosceles triangle 



exact 

Polya (|y_ 
Freitas (pL5l 
Protter ( 
Theorem 
Theorem 



30.7054 
26.3189 
23.0291 
19.0338 
23.9381 
29.8449 



(b) Half of the equilateral trian- 
gle 



Table 1: Triangles with known first eigenvalues 





arm 2 


arm 4 


Polya (1. 
Freitas (1 
Protter ( 
Theorem 
Theorem 
Upper bo 


11 

.15 

1.14 

1.3 

1.4 

unc 


I (6.1) in fisl 


23.5404 
20.3972 
17.0662 
20.9906 
27.0781 
27.6695 


11.4865 
12.4937 
12.8437 
12.9675 
18.8754 
18.9749 



Table 2: Tall isosceles triangles 



values of the lower bounds for an isosceles triangle with base 1.95 and arms 1 are 
given in Table [3] 



Polya (1. 
Freitas (1 
Protter ( 
Theorem 
Theorem 


11 
.15 

i.U 

1.3 
1.4 


1 


105.206 
210.273 
205.698 
212.735 
185.161 


Conjecture 1.2 in 18 




Lower bound 
Upper bound 


251.077 
299.7 



Table 3: Wide isosceles triangle 



Clearly Freitas, Protter and Theorem 1.3 
values are not very accurate 



are the best in this case, but those 



In fact, neither bound allows us to prove the second 
part of the Theorem |1.1| in the case of obtuse triangles. Note that Conjecture 1.2 
18 could provide a lower bound strong enough for this task. Theorems 2 and 



5.1 in 11 could also provide good enough bounds. 
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3. Symmetrization techniques 

In this section we present various geometric transformations that decrease the 
first eigenvalue. The following subsections contain strict definitions of three dif- 
ferent kinds of symmetrization. Here we just remark that the important common 
property of those transformations is that they are contractions on T4^^'^(M'')nCc(M'') 
and isometrics on L^(]R'*). Those properties, together with the minimax formula 
for the first Dirichlet eigenvalue, give Ai(ri) > Ai(ri*). The most general reference 
for those results is [12) . 



3.1. Steiner Symmetrization. We start with the well known Steiner symmetriza- 
tion (see for example |15j Note A] or 12 Chapter 2]). Fix a line I = {ax+b : a; € M} 
where a and b are arbitrary points on the plane. Let {lt}teM. be a family of lines 
perpendicular to I such that for each t the line It passes through the point at + b. 
For an arbitrary domain f2 we define its Steiner symmetrization f2* with respect to 
I as a convex domain symmetric with respect to / and such that for every line It 



(3.1) 



\n*nit\ = \nnitl 



where | • | denotes the 1-dimensional Lebesgue measure. 

Intuitively, we look at the cross-sections of perpendicular to I and we center 
them around / (see Figure |3] for an example) . 

This procedure has many interesting properties (see 12 15 ). While the area 
remains fixed, the perimeter decreases and the inradius increases. But the most 
interesting property from our point of view is that the first eigenvalue of the Dirich- 
let Laplacian on f2 is bigger than that on fl*. This basic property has been widely 
used in the proofs of many isoperimetric bounds for eigenvalues in many settings, 
see for example [ll[l2 15 . We will use this property to prove our results. 



3.2. Continuous Steiner Symmetrization. The second type of symmetrization 
is the continuous Steiner symmetrization introduced by Polya and Szego in [l5|Note 
B] . Different versions of it have been studied by Solynin '20] and Brock [4j[5j. The 
difference in case of convex domains is only in the way of defining "continuity 
parameters" . The version studied in [4jj5] is more general and it works even for not 



connected domains. We refer the reader to these three papers and to 12 for the 
properties of this transformation. Here we only give the definition valid for convex 
domains. As above we look at all the intervals (at, bt) which are the intersections of 
with /(. Let {a't, b't) be the Steiner symmetrized interval. Let af = at +a{at — at) 
and bf = bt + aib't — bt). Hence we are shifting the intersections with constant 
speed from their initial position to the fully symmetrized position. We define the 
continuous Steiner symmetrization f2" of a domain by 

(3.2) n'^ = \J{atX)- 

Figure [3] shows the action of continuous Steiner symmetrization on triangles. 
We can see that = Q and fi^ = fi*. The first eigenvalue of the domain fi" is 
decreasing when a is increasing (see [4j|5][l2]) just like in the case of the classical 
Steiner symmetrization. Note that our "continuity parameter" a is related to the 



one in 12 by 1 — a = e 



-t 
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Figure 3: Continuous Steiner symmetrization 



3.3. Polarization. The last technique we want to introduce is called polarization. 
It was used in |9 to prove general inequalities for capacities and eigenvalues. It has 
also been useful in studying other types of symmetrization (see [6][20]) and heat 
kernels for certain operators (see |8]). Let I be an arbitrary line, Hi and H2 two 
halfspaces with boundary I. For x £ H2 let x denote the reflection of x with respect 
to I. Polarization of a domain is defined pointwise by the following transformation 

Definition 3.1. The polarization of a set f2 with respect to I is a set with the 
following properties. 

(1) Ifx ennHi then x € fl^ . 

(2) If x^nn H2 then x e fl'" . 

(3) // both X and x are in ^l, then both x and x are in il^ , 

In essence, we reflect the set H, with respect to /, the intersection of H. with its 
reflection is in the polarized domain H,^ and all other points of Hi that are in either 

or its reflection also belong to fi^. Figure [4] shows a triangle, its reflection with 
respect to I and the polarized domain (outlined polygon) . 

This simple procedure has many interesting monotonicity properties, among 
them is the monotonicity of the eigenvalues. For this, and many other proper- 
ties we refer the reader to l6|p] and to references in these papers. 



4. The proofs of the lower bounds 



Before we prove Theorems 1.3 and 1.4 let us present a one simple application 
of the continuous Steiner symmetrization. This and other kinds of symmetrization 
were introduced in Section |3] 

Let T be an acute isosceles triangle with area A and arms of length d. Let I be 
a line perpendicular to one of the arms and passing through the midpoint of this 
arm. If we apply the continuous Steiner symmetrization with respect to I we get 

Lemma 4.1. Amongst triangles T with diameter d and area A, isosceles trian- 
gle with arms d maximizes the first eigenvalue and isosceles triangle with base d 
minimizes it (See FigurelSy. 
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Figure 4: Polarization of sets 



This simple result shows that upper bounds for the eigenvalues can also be ob- 
tained using symmetrization techniques. Using the largest circular sector contained 
in T we can get very accurate upper bounds for the eigenvalues of triangles. 



We can use polarization to give an alternate proof of (1.151. A very similar 



procedure, although more complicated, will be used later to prove Theorem |1.4[ 
Let T be an arbitrary triangle with the longest side of length d. First, we use the 
Steiner symmetrization with respect to a line perpendicular to the longest side. As 
a result we have an isosceles triangle ABC with the base of length d. 




Figure 5: Freitas's bound 



The further construction is shown on Figure [5) We divide the longest side AC 
into four equal parts. Then we construct the line EH perpendicular to AC . Finally 
we construct the bisector li of the angle AGH. We apply the polarization with 
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respect to the line li. The darker triangle with vertex H is the result of the 
reflection with respect to ^i. Note that we need the line li to cut the side AC 
between the points E and F. If not, then we would not be able to symmetrize the 
other half of the triangle (the bisector would cut the sector EG, giving unnecessary 



reflections). Fortunately one can check that under the condition (2.7 1 this is always 
true. 

Now we perform one more polarization with respect to the line I2, which is the 
bisector of the angle HIB. As can be seen on Figure [5] the left side of the triangle 
has been changed into a rectangle. 

If we repeat the procedure on the other side we get the rectangle with the base 
EE and the height BD. But for a rectangle we have the following explicit formula 
for the first eigenvalue. 



(4.1) 



^l\Rectangle = TT^ (fl ^ + 5 ^) 



where a and h are the lengths of the sides. In our case the lengths are \EF\ — d/2 
and \BD\ = h. This gives Freitas's bound (1.15). 



To get a sharper bound we need to symmetrize the triangle into a rectangle with 
a longer base and a shorter height. We start with the same Steiner symmetrization 
as before, That is, we symmetrize with respect to the line perpendicular to the 
longest side. Next, we perform one more Steiner symmetrization but with respect 
to the longest side. This gives a rhombus with diagonals of length \CD\ — h and 
\EF\ = h (see Figure |6]). 




Figure 6: Theorem 1.3 



If we apply the Steiner symmetrization one more time but with respect to the 
height DE of the rhombus, we obtain the rectangle with base AC and height DE. 
Using Pythagorean theorem we find that \AC\ = ^/cP + h^ /2. Since the area A of 
the triangle remains constant under symmetrization, we also have \DE\ = A/\AC\. 
These, together with (4.1 1, give the proof of Theorem 1.3 Note, that the proof of 



Theorem 1.3 is actually easier than the proof of Freitas's bound (1.151 



It remains to prove Theorem |1.4[ Let us begin with a lemma. 

Lemma 4.2. Let T be a triangle with the smallest angle 7. Let T' he a triangle 
with same area A and same smallest angle 7 hut with a smaller diameter. Then 



(4.2) 



AiIt > All 
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Remark 4.3. What this lemma essentially says is that we can continuously deform 
a triangle into an isosceles triangle, preserving the smallest angle and the area in 
such a way that the first eigenvalue will be decreasing. This immediately gives the 
first inequality in Theorem\1.4\ 



Proof. This is perhaps the most compHcated appHcation of any symmetrization 
techniques. We have to apply a suitable sequence of polarizations to obtain the 
result. First, let be a triangle similar to T' but with area A + e. We will 
symmetrize T into a set contained in T^. Then 

(4.3) Al|T>Al|T^ 

But when £ ^ 0, the eigenvalue of converges to the eigenvalue of T' and this 
ends the proof. Indeed, the triangles are similar hence we get the convergence due 
to scaling property of the eigenvalues. 

Ideally, we would like to just define a sequence polarizations that folds T inside 
T^. Unfortunately, it is virtually impossible to choose a correct line even for the 
first polarization. Each polarization must be with respect to a line cutting the 
longest and the shortest sides, but its exact position is not clear. The first line 
should be close to the vertex and the following lines should move away from the 
vertex. To give a precise position of each line we need to consider a temporary 
reversed sequence of transformations (very similar to the sequence of polarizations 



performed in the proof of 1.151. The first line in this reversed sequence (or the last 
line for polarization) is easy to define. Having this line we can get another, and so 
on. The precise construction is split into 4 steps. 

Step 1: Definition of elementary transformation. 

We start with two triangles T — ABC and — DBE such that E is on the 
interval EC and A is on the interval DB. We also assume that the angle with 
vertex B is the smallest in the triangle T and that the area of T'^ equals the area of 



T plus e. The triangles are shown on Figure 7a Let li be the bisector of the angle 



between the intervals DE and AC . If we reflect T with respect to this bisector we 



obtain another triangle A' B'C shown on Figure 7b Let G be the intersection of 
C'B' and AC. 

We define the elementary transformation of T to be the polygon ABEC'G. This 
construction is valid due to two conditions. The first, the area of DBE is bigger than 
the area of ABC . The second, the angle ZACB is smaller than the angle ZEDB. 
Those two conditions guarantee that the bisector cuts the sides AC and BC and 
not the side AB and validates the construction of the elementary transformation. 

The same transformation can be applied to more general domains. For example 
any domain contained in an infinite cone ACB with vertex at C. We will need this 
observation later on. 

Notice also that the elementary transformation is equivalent to polarization with 
respect to the bisector for any domain containing ABC as long as: the intersection 
of this domain with the bisector is the same as the intersection of ABC with the 
bisector, the new part (triangle C'GF in the case above) does not intersect the 
domain. In particular we could take a domain consisting of ABC and any subset 
of a half-plane disjoint with ABC and with boundary containing AB unless this 
domain intersects the bisector or C'GF. We also do not need C'GF to be a triangle, 
it is enough that this set does not intersect the domain. 
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Figure 7: Triangles T = ABC and = DBE. 



Step 2: Sequence of elementary transformations. 

Now we want to describe a sequence of elementary transformations that changes 
T = ABC into a subset of = DBE. 

The first transformation is described in Step 1. It is clear that the area of ADF 
is bigger than the area of C'CF. We have 

(4.4) ZGC'F = /LACB < ZACB + ZABC = ZDAF. 

Therefore the second elementary transformation can be applied to triangles C'CF 



and ADF (see Figure 8a I . Here we completely disregard the presence of the quadri- 
lateral AFEB. 

This elementary transformation introduces a new triangular piece that may inter- 



sect the triangle ABC near vertex A (see Figure 8b I . The intersection means that 
the elementary transformation is not equivalent to the polarization of ABEC'G 
(we did not care about the quadrilateral AFEB while making a transformation). 
This intersection will not occur if the triangle ADF is large enough to fit the tri- 
angle C'CF inside. It would follow that the second elementary transformation is 
also equivalent to polarization. In such a case the construction of the sequence of 
transformations would be finished. 

If the intersection occurs, we apply another elementary transformation to the tri- 
angle introduced during the previous transformation and to the difference between 
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(a) bisector (b) transformation and intersection 



Figure 8: The second transformation 



the triangles considered for the previous transformation (triangle AGH). Here we 
again disregard the presence of other part of the domain. The third elementary 
transformation is shown on Figure |9] To validate this step we need to check that 
the angle ZAGH is bigger than the angle ZACB. Indeed, it is equal to the angle 
ZGC'F{= ZACB) plus the angle ZC'FG. 




Figure 9: The third transformation 



If the newly introduced triangular piece does not intersect the already obtained 
domain (as on Figure [9]), the construction in finished. If it does, we apply another 
elementary transformation (or as many elementary transformations as needed) to 
fully "fold" the triangle AEG inside of the triangle DEE. The final subset of DEE 
that has the same area as AEG will have a spiral-like pattern of triangular pieces 
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introduced by consecutive elementary transformations. First we need to show that 
this process is finite. Later we need to modify the procedure to get a sequence of 
vahd polarizations. In particular we need to avoid self-intersecting that happens at 
every stage. 

Step 3: Finiteness of the construction. 

To estimate the area of the subset of ADF that is covered using three consecutive 
elementary transformations, we shift the triangular pieces obtained each time so 
that the angle equal to the angle ACB has a vertex at D, A and F respectively (after 
the second, third and fourth elementary transformation). The triangle A'B'C' is 
shifted so that vertex C" moves to D, G moves to F' and the angle ZFDF' is equal 
to the angle ZACB. We do the same with other triangles. See Figure [TO] for the 
picture of the shifted pieces. This decreases the covered area, but does not change 
the angles between any of the lines. 




Figure 10: Covering using the angle ACB (arrow) 



We have 

(4.5) ZA'F'D' + ZAF'A' = ZF'DF + ZF'FD. 

This implies that the angle ZA'F'D' equals the angle ZAFD. Similarly we show 
that other angles of the triangles AFD and A'F'D' are also equal. This means that 
three elementary transformations reduce the uncovered part to a triangle similar to 
AFD. The size of this triangle is not bigger than the size of A'D'F'. Therefore the 
area of the uncovered part is shrinking at a constant rate. Due to the difference in 
the area of and T the procedure must end as soon as the area of the uncovered 
part is smaller than e. 

Step 4: Sequence of polarizations. 

The presence of intersections in the second and all consecutive elementary trans- 
formations stops us from using the sequence of elementary transformations as a 
sequence of polarizations. To fix this problem we consider the same sequence of 
transformations, but in the reverse order. Each elementary transformation was ap- 
plied to a triangle that is a part of the reflection of the triangle ABC. Therefore 
we can treat those elementary transformations as a transformations on ABC. 



16 



BARTLOMIEJ SIUDEJA 



As we remarked in Step 1 such transformations are equivalent to polarizations. 
We use the reflection lines used by the elementary transformations to produce po- 
larizations. Figures | lla| and | llb| show the domain after the first two polarizations 
or the last two elementary transformations (outlined sets). The construction from 
Step 2 ensures that the triangular pieces introduced at every step have no inter- 
section with the reflection lines. We also avoid self-intersecting since the triangular 
pieces are introduced in the reverse order. This means that we have a sequence of 
valid polarizations. 

The last polarization (the first elementary transformation) fits the whole domain 
inside T'^ = DBE (see Figure 12). If the sequence is longer than on the example, 
we proceed in the same manner, building the spiral-like structure starting from the 
smallest inner triangle. 

This proves that for arbitrary e > there exists a finite sequence of polarizations, 
which transforms T into a subset of T^. 

This sequence of polarizations can be constructed even for domains consisting 
of triangle ABC and any set contained in the half-plane with boundary AB and 
disjoint with ABC . 





(a) The first polarization (tlie third (b) The second polarization (the second transfor- 
transformation) mation) 



Figure 11: Reversed steps 



□ 



Remark 4.4. The same proof as for Lemma \4-2\ works for any domain contained 
in an infinite cone ACB and containing T . 

Lemma 4.5. Let CI be a kite, that is a quadrilateral symmetric with respect to the 
longer diagonal and with perpendicular diagonals. Assume also that the length of 
the longer diagonal is smaller than the length of one of the sides. Given the fixed 
area and the smallest angle, the first eigenvalue is decreasing with the diameter. 
Ln particular, an isosceles triangle has a maximal first eigenvalue among kites with 
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Figure 12: The third polarization (the first transformation) 



fixed area and the smallest angle, 
minimal first eigenvalue. 



the kite consisting of two isosceles triangles has 



Proof. We can split an isosceles triangle ABC into two right triangles and use 
Lemma |4.2| and the remark above to symmetrize any one of those right triangles 
(for example DBC on Figure 13a). Furthermore, we can perform all but the last 



polarization on both right triangles (Figure 13b shows the domain at this stage of 
construction) . Finally we can perform the last polarizations to get the quadrilateral 
with vertex E. The last transformations are indeed polarizations since the bisectors 
are cutting the line BD below the point B (certainly not between points D and 
E). 

□ 



The proof of the second inequality in Theorem 1.4 requires a repeated use of 
the above results. As before we take a slightly larger sector S{A + e,^) instead of 
S{A, 7). We need to transform the isosceles triangle I{A^ 7) into a subset of S{A + 
£,7). Using Lemma 4.5 we can symmetrize the triangle ABC into a quadrilateral 



(see Figure 14a) with longer diagonal equal to the radius of the circular sector 
5(A + £,7). 

Now we divide both halves of the quadrilateral into N (to be chosen later) parts 
by dividing the angle with vertex B into equal parts (see Figure 14b[ ). 

We work on each half separately. Consider a triangle formed by all but the inner 
most part. We use Remark |4.4| to symmetrize this triangle into a "more isosceles" 
triangle with one vertex on an intersection of a circular part and the innermost 
line. This makes the part outside of the circular sector smaller. Now we repeat this 
using decreased number of parts. The picture after two steps is shown on Figure 
[TS) By choosing a large enough N and performing all steps we can fit the whole 
half of the quadrilateral inside of the circular sector. 
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(a) one-sided symmetrization (b) two-sided symmetrization be- 

fore the last step 



Figure 13: Isosceles triangle 




(a) quadrilateral (b) subdivisions 

Figure 14: Symmetrized isosceles triangle 



This ends the proof of the second inequality in the Theorem 1.4 The last thing 
to prove is the monotonicity property for the isosceles triangles. Here we need to 
apply the continuous symmetrization to obtain the result. Suppose that we start 
with the isosceles triangle I{A, 7) with the vertex angle smaller than tt/3. We want 
to show that if we increase the angle to 7r/3 > 7' > 7 while the area remains 
fixed, then the first eigenvalue decreases. First, we apply the continuous Steiner 
symmetrization with respect to the line perpendicular to one of the arms. The 
length of the base increases, and its maximal length is obtained when we reach the 
full Steiner symmetrization. Compare with Figure |3] where the shortest side of fi" 
increases with a. 
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If this maximum is smaller than the base of /(A, 7'), we use the Steiner sym- 
metrization to get an isosceles triangle with the angle 7" between 7 and 7' and the 
base equal to this maximum. We can repeat the above procedure on the triangle 
7(A,7"). 

On the other hand, if the maximum is bigger than the base of J(A, 7'), then 
we can stop the continuous Steiner symmetrization at the time when the length of 
the enlarged base of I{A, 7) is equal to the base of I{A, 7') (on Figure |3] we choose 
a such that the shortest side is equal to the base of 7(^,7'). Now, the Steiner 
symmetrization with respect to the base gives an isosceles triangle with the base 
equal to the base of 7(^,7'). This implies that it must be equal to 7(^,7'). 

Suppose that the vertex angle 7 > tt/3. The same procedure as above shortens 
the base in this case. The same argument applies, but with maximum replaced 
with minimum. This ends the proof of the Theorem |1.4[ 



5. The proofs of the upper bounds 

In this section we prove Theorem The approach we take is based on vari- 
ational bounds for the second eigenvalue. This approach is a modification of the 
similar method used in [18 . We start with the minimax formula for the second 
eigenvalue. 



(5.1) A2ID = inf sup 



|V(/i + a/2 



12 



12 



where /i,/2 G Hq(D) are linearly independent. This formula is a special case of 
the general minimax formula for an arbitrary eigenvalue (see e.g. [?]). As in 18 



we will use known eigenfunctions for equilateral or right triangles to obtain test 
functions for arbitrary triangles. 

Consider the equilateral triangle Tf. with vertices (0,0), (1,0) and (1/2, -s/S/Z). 
The complete set of eigenfunctions is well known. For the exact formulas for these 



eigenfunctions we refer the reader to 13 , 14 . In particular, [13] gives the simple 
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formulas that we will use in this paper. Let 

(5.2) ,= l(^2x-l), 

(5.3) t = 7r(l-^y 



V3- 

The first eigenfunction is given by 

(5.4) y) = (cos(3z) — cos(i)) sin(t). 

The second eigenvalue has multiplicity two. We will follow the notation from 



13 to name these eigenfunctions. All eigenfunctions can be divided into two kinds: 
symmetric and antisymmetric. The first will be denoted by S followed by two num- 
bers identifying the eigenfunction. The second by A also followed by two numbers. 



For the details about this notation we refer the reader to 13 . The two eigenfunc- 
tions belonging to the second eigenvalue are 

(5.5) '■PS2i{x) = cos(4z) sin(2f) -I- cos(5z) sin(i) — cos(2;) sin(3t), 

(5.6) '4^A2i{x) = sin(4z) sin(2t) + sin(5z) sin(t) — sin(z) sin(3t). 

Let T be an arbitrary triangle. We can assume that one side of this triangle is 
equal to the segment from (0,0) to (1,0), and that the last vertex (w,f) is in the 
upper halfspace. Then, there exists a unique linear transformation L from T onto 
Tg. As in |18j, we will compose L with the eigenfunctions of Tg to obtain suitable 
test functions for T. 



Using formula (5.11 we obtain an upper bound 

(5.7) A2|T<sup^y-^ 

where ga is a linear combination of two known eigenfunctions. 
As the first two test functions we can take 

(5.8) gi{x, y) = (ps2i{x, y) + aipsii{x, y), 

(5.9) 9l{x,y) = (pA2i{x,y) + a(psii{x,y). 

If the triangle T is almost equilateral, we can expect that its eigenvalues and 
eigenfunctions are similar to the eigenvalues and eigenfunctions of T^. Also, the 



linear transformation L should not perturb the bound (5.7) significantly. Therefore, 
we can expect that for almost equilateral triangles this should be a good upper 
bound for A2|t- 



Notice that the linear combinations in (5.8 1 and (5.9 1 consist of two orthogonal 
functions. Hence the second norm of this combination is just the sum of the second 
norms. Therefore 



(5.10) sup 



/j,|V(5„oL)|2 _ aa^ + ba 
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where a,c,e,f are strictly positive. This rational function has a limit a/e, as a 
±cxD. By taking its derivative we can also find the critical points 



(5.11) 
(5.12) 



ai = 



1 

bd 



if-ce^ l^4b^ef+{2af-2ce)A , 



a2 



^ [af - ce + ^ v/462e/ + (2a/ - 2ce)' 



If we evaluate the function at those points and simplify, we get 



(5.13) 
(5.14) 



Ci = 



Co 



1 

1 

2^f 



af + ce^\ ^Ab^ef + (2a/ - 2cey 



af + ce+^^4b^ef- 



(2a/ - 2cey 



The expression under the root is clearly nonnegative. It is zero if 6 = and af — ce. 
In such case Ci = C2 = a/e, hence the maximum is a/e. If the expression under the 
root is positive, then we have two distinct critical points, and C2 > Ci. This means 
that this rational function has an absolute maximum C2 and absolute minimum 
Ci. 

This leads to a new bound for the eigenvalue 



(5.15) 



A 



2\T 



1 

2^f 



af + ce+^ ^Ab^ef + (2a/ - 2ce)' 



To finish the proof of Theorem |1.1| we need to use one of the lower bounds for 
the first eigenvalue. Due to the simple form of the bound, we will use Freitas's 
result ( |1.15 1 . To prove the first bound in Theorem |1.1| we need to show that 

167r2 



(5.16) 



((5.15) - (1.15))i?^ < 



27 



where we have used the reference number of the equation for the value given by it. 
The number on the right side of the inequality is the exact value obtained for the 
equilateral triangle. 

The second bound will be proved if we can show that 



(5.17) 



< 



7 



Notice that this time we use Polya 's isop erimetric bound (1.1 1. 

Let us begin with the proof of (5.161. We first use 5^ as a test function. The 
expressions (5.15 1 and (1.15 1 can be written in terms of the vertices of the triangle 
T. But we assumed that the vertices are (0,0), (1,0) and {u,v). We denote the 
lengths of the sides of the triangle by 1 , M = + w^, N — — uY + f ^ and 
we can assume that iV > M > 1. Then the bound ( 5.16| ) (in terms of the lengths 
of the sides) is equivalent to 

> -1612800A^2(1 + M + NfTi"^ + 27|-413343iVV 
(5.18) +11200(9(Af^ - 1)^ + 2{AP + l)N^ + 2QN^)tj:^ 

+iV2 \/655128046899y2 _ 740710656001/M^7r2 + 80281600001^27^41 ^ 

where V + -2 and W = + + 1. 
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The expression is quite complicated since the integrals of the function gl^ are 
very cumbersome to calculate. This task could be accomplished by hand since this 
function is just a sum of the products of the trigonometric functions. However, 
symbolic calculations in Mathematica are used to obtain this expression in a short 
time. The script performing all the calculations from this section is included in the 
last section. 

One can check numerically where this inequality is true. If we put U = N — M 
then the set of all possible triangles can be characterized by G [0, 1) and M > 1. 



The dark part of Figure 16 corresponds to those triangles for which the inequality 
is valid. 




Figure 16: Bound 1/Case 1: The numerical solution 



It is clear that not all triangles can be handled this way. Hence we have to use 
more then one test function, and divide all triangles into subregions of U G [0, 1) 
and M > 1. We will define other test functions at the end of the section. 



We now prove inequality (5.181 on the gray rectangle shown in Figure 16 More 



precisely we take U € [0, 0.03] and M € [1.03, 1.39]. 



The inequality (5.181 can be written as 

(5.19) P{N, M) + Q{N, M) y/R{N, M) < 0, 

where P, Q and R are polynomials in N and M. It will be proved, if we can show 
that 

(5.20) P{N,M)<0, 

(5.21) Q^{N, M)R{N, M) - P'^{M, N) < 0. 

This is a system of polynomial inequalities. Unfortunately the degrees of those 
polynomials are 4 and 8 in each variable. Therefore there is almost no hope to 
solve this system using any conventional method. 

Instead, we developed an algorithm for solving such polynomial inequalities on 
rectangles. The next section contains a detailed description of this algorithm and 
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the proof of its correctness. It turned out that in our case this method gives the 
proof of the inequality for any test function we tried. 

To finish the proof of the first part of Theorem |1.1| we need to define the other 
test functions and rectangles for each corresponding inequality. 



We already have two test functions (5.8) and i5.9\ that came from the equilateral 



triangle. We can also use the first two eigenfunctions of the half of the equilateral 
triangle. That is, of the right triangle with angle 7r/6. Its eigenfunctions are certain 
antisymmetric eigenfunctions of the equilateral triangle. Hence, we get the third 
test function 

(5.22) gl{x, y) = (pA3i{x, y) + a(pA2i{x, y), 



where (using notation (5.2 1 and (|5.3[)) 



(5.23) fA3i {x,y) — sin 5z sin 3t — sin 2z sin At — sin 7z sin t. 

In this case we also need to get a new linear transformation L' that transforms the 
given triangle into the right triangle. 

The last case when all eigenfunctions are known is the right isosceles triangle. 
Since this triangle is "half" of a square, its eigenfunctions are equal to eigenfunctions 
of the square with diagonal nodal lines. We get 

(5.24) gl^{x,y) = (l)2{x,y) + a(j)i{x,y), 
where 

(5.25) (j)i{x,y) = sin27ra;sin7ry + sin7ra;sin27r?/, 

(5.26) 4'2{x,y) = sin37rxsin7ry — sin7ra::sin37ry. 

We can also mix the eigenfunctions from the different triangles provided we use 
appropriate linear transformation for each of them. In this manner we obtain the 
last (fifth) case needed to prove the theorem by taking gl^{x,y) and ^ip2{x,y). 

Each of the five cases requires a rectangle on which we can prove the bound. We 
take 

(1) U e [0,0.03], M € [1.03,1.39], 

(2) U e [0,0.2], Me [1,1.03], 

(3) U e [0,1), M € [1.39, oo), 

(4) Ue [0.2,1), Me [1,1.39], 

(5) U e [0.03,0.2], Me [1.03,1.39]. 

Notice that these rectangles exactly cover the infinite strip [0, 1) x [1, oo). Since the 
sides of the triangle are A'' > M > 1 and U = N — M, the strip includes all possible 
combinations of lengths of the sides. 

The proof in cases (3)-(5) is exactly the same as in case (1). Additional step 
has to be performed in case (2). Consider a triangle T' similar to T with sides 1, 



N' = M/N, M' = 1/N. We have 1 > N' > M'. Consider inequality (5.161 for 



T' . Note that the diameter of T' used in (1.15) is now 1 {N in other cases). Just 



like before we get an expression similar to (5.18) but involving N' and A/'. After a 



change of variable M' —>■ N' — M' + 1 and N' ^ I — U' we can apply our algorithm 



with a rectangle given in case (2). This proves the bound (5.16) for T', and hence 
for T since the bound is invariant under scaling. 
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We need to check that the transformation 
(5-27) — ^ = l-U', 

(5.28) — ^ — - M', 

^ ' M+U 

changes the rectangle U € [0, 0.2), Me [1, 1.03) into a subset of the same rectangle 
in primed variables. From the first equation we get 

(5.29) 

From the second equation 

This ends the proof of the first part of Theorem The second part requires an 
additional argument. We still need to use test functions g\{x,y) through g^{x,y) 
on the following rectangles 

(1) U e [0,0.09], M e [1,1.37]; 

(2) not needed; 

(3) U e [0,0.42], M e [1.37,2.05]; 

(4) U e [0.09,0.2], Mg [1,1.37]; 

(5) U e [0.2,0.42], M € [1,1.37]. 

Notice that in each case U < 0.42. Since the triangles are acute we have iV^ < 
+ 1, hence U < \/2 — 1. Also, all these rectangles cover only the cases with 
M < 2.05. 

Therefore we need a circular sector type-bound similar to the one in |18| . We 
use the sector bound from Theorem 1.4 for the first eigenvalue and the upper bound 
for the second eigenvalue based on the biggest circular sector contained in the given 
triangle. Suppose that we have an acute triangle with area A and smallest angle 7. 
Let N > M > 2.05 be the two longest sides. Since the triangle is acute, the largest 
sector contained in the triangle has the radius equal to the longest altitude H. We 
get 



(5.31) ^ 



Triangle 

This bound is certainly good enough, but it is hard to deal with due to the 
complicated formula for the radius of the circular sector S{A,'-f). Therefore we 
have to rely on the weaker bound 

A2|s(7_H2/2,7) 



(5.32) ^ 



< 



Triangle ^l\s{'rMN/2,'r) 

This bound also follows from Theorem |1.4| by taking the smallest sector containing 



the isosceles triangle /(A, 7) (first inequality in Theorem 1.4 1 with the arms of the 
length y/MN. 

The eigenvalues of the circular sectors are given in terms of the zeros j^^n of the 
Bessel function of index v. Here n indicates n-th smallest zero. We have 



(5.33) Ai|s(7fl2/2,7) = ^ J7r/7,1' 



(5.34) A2|5(7_f/V2,7) ^ •^i2('5'ector) = i? 



7r/7,2. 
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Notice that we only have the inequahty for the second eigenvalue due to the pres- 
ence of another eigenfunction X2i{Sector). This eigenvalue may be smaller than 
Xi2{Sector). 

We can use the bounds for the zeros of the Bessel functions proved in |17j. That 

is, 

(5.35) ^ - < < " - M^'H' 

where Uk are zeros of the Airy function with ai ~ —2.3381 and 02 ~ —4.0879. 
Using the last two facts we obtain 



A2 
Ai 

(5.36) 



< 



NM I ^'-fj^+lA"^ 



,3/7T \ 2 



X 



Triangle \ \/2 

_NM 



2 



where x — (71/7)^/^. Note that for acute triangles we have N"^ < We need to 

consider two cases based on the length of H. First assume that > — 1/16. 
This is equivalent to the condition that the altitude H divides the shortest side 
(with the length 1) into two parts with one of them not longer than 1/4. Using the 
inequality < + I, we obtain 

NM M^/AP + 1 I — -( 1 

If AP - 1/4 < i?^ < - 1/16 then the altitude H divides the shortest side 
into two parts, one of length (5 satisfying 1/4 < 8 < 1/2 and the other of length 
1 - 6. Then 

(5.38) N^^H^ + {1- = AP - (5^ + (1 - 5f = AP + 1- 2S < AP + 
This gives 

Let z{M) = max{zi(M),^2(-M)}- Note that z(M) is decreasing since both zi(M) 
and Z2{M) are decreasing. Making the substitution y = x~'^ we get 



(5.40) 



A2 
Ai 



< z{M) 

Triangle 



2 2^ ^ 
20 "'2i/ 
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where the constants Ci satisfy 

(5.41) = 0.83133 >0. 

(5.42) c..-jj-. 1,70183 >0, 



lOaf - 10aia2 + 3a| 



(5.43) C3 = — 1;^^-2 ^ o^^ggg^ ^ 0^ 



(5.44) C4 = --^ w 1.85575 > 0. 

v2 



We want to show that the right hand side of (5.401 is decreasing with M (note 
that y also depends on M since it depends on 7). First we can show that this 
expression is increasing with y for a fixed Af . Indeed, the derivative with respect 
to 2/ is C2 — (i4-c4y)^ positive for y > Q due to the condition C2 > C3. This 

means that the expression is also increasing in 7 since y = (71/7)^^/^. If we fix M, 
then from all triangles with sides 1, M > 1, > M , the isosceles triangle (Af = N) 
has the biggest angle 7. In this case cos 7 = 1 — 2W^- ^ result we get an upper 
bound for the ratio of the first two eigenvalues in terms of M. That is, we have 

2 



(5.45) ^ 



< z{M) ( ci + C2yM + 



Triangle \ ^ + ^^^VM , 

with = (arccos(l — (2Af^)~^)/7r)^/'^. But yM and z{M) are decreasing with M, 
hence the right hand side is decreasing with M . To finish the proof we just need 
to check that this is smaller than 7/3, as required in ( |5.17 ), for M = 2.05 (we get 
w 2.3285 < 7/3). This ends the proof of Theorem [TT] 

6. An algorithm for polynomial inequalities 

This section gives the algorithm for proving polynomial inequalities in two vari- 
ables with arbitrary degrees. The domain we deal with is a rectangle. The proof of 
the correctness is also given. 

Suppose that we have an inequality 

n m 

(6.1) P(a:,2;) = ^^Q^,xV <0, 

for X E (0, a) and y E (0, b). Any other rectangle can be shifted to the origin, hence 
reducing the problem to this case. 

The idea behind the algorithm is very simple. For any monomial 

(6.2) c,jx'y^ < c,j mm{ax'-'^y^ ,bx'y'''^}, if c,j > 0, 

(6.3) Cijx'y^ < c,^jma.x{a-^x'+^y\b-^x'y'+^}, if j < 0. 

We can use this simple observation to reduce the number of positive coefficients in 
P{x,y). Clearly, if we apply any of the above inequalities finite number of times on 
any of the monomials in P(x,y), we obtain an upper bound for P(x, y). If we can 



reduce the whole polynomial to 0, we proved inequality (6.1 1. 

We need to describe a sequence of reductions, that leads to a constant or to a 
polynomial with only positive coefficients. If we get at least one strictly positive 
coefficient we cannot say that inequality is false, since we are using an upper bound 
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for P{x,y). Thus, the algorithm does not always work although it is possible to 
generalize it. In case of the false answer we can divide the rectangle into four 
identical sub-rectangles, and rerun the algorithm on each of them. As long as the 
inequality P{x,y) < is strict, this recursive procedure should give the proof. 

It is worth noting that for 8 out of the 9 polynomials we have in the previous 
section this method works on the whole rectangle (without sub-dividing). In the 
case of the third test function in the second part of Theorem |1.1| we need to split 
the given rectangle into halves. The gray rectangle on Figure [16] in the previous 
section is almost as big as possible given that the inequality is true only for the 
dark points. Hence the method works very well in this case. By running the script 
from the last section one can see that in all other cases the method is also very 
efficient. 

Now we will define the optimal sequence of reductions. The only way to reduce 
a positive coefficient is by lowering the power of one of the variables (using ( |6.2[ )). 
Similarly, to reduce a negative coefficient we have to increase one of the powers 



(using (6.3l). Each time, two of the coefficients combine giving a new, possibly 



negative coefficient. Write 

n 

(6.4) P(x,y) = ^a;^Q,(2/), 

4=0 

where 

(6.5) Qj(y) = ^Qj-y-'. 

3=0 

To avoid ambiguity, we start with Q„. Any negative coefficient in Qn can be 
used only to combine with some positive coefficient with higher power oi y. It is 
impossible to use them to interact with Qi for i < n. Therefore we inductively 
(starting from j = 0) check if c„j- < and in case this is true we get an upper 
bound 

(6.6) c„,,x"j/^' + c„.j+ia;"y^+i < (c„.,6-i + c„j+i)x"y^'+\ 

and we redefine c„j+i = c^ jb^^ + c„ j+i. If the last coefficient turns out to be 
negative, we can just change it to 0. 

As a result we changed Qn{y) into a polynomial Qniu) with nonnegative co- 
efficients (possibly all equal to 0). A positive coefficient can only be altered by 
lowering one of the powers. We could lower the power of y but, ultimately, if we 
want to obtain a constant as a final bound we have to also lower the power of x in 
all coefficients of x'^Q'^{y). Hence we get an upper bound 

(6.7) x^Q[,{y) + x"'iQ„_i(y) < {aQ'^{y) + Q„_i(y))x"-\ 

and we redefine Qn-i{y) — o,Q'n{v) + Qn-iiu)- This approach guarantees that any 
negative coefficient of Qn-i{y) can be used to reduce as many positive coefficients 
from Qn{y) as possible. 

Now we repeat the whole procedure for Qn-i{y) and for all others by induction. 
At the end we get (9o(y) which has only nonnegative coefficients. If any of those 
coefficients is strictly positive, the algorithm failed. But Qo{y) = means that 
P{x,y) < on the rectangle (0,a) x (0,6). 
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The implementation of this algorithm is a paxt of the script in the following 
section. The following simple example shows how the algorithm works. Let 

(6.8) P{x, y) = x^y^ - x^y + 2xy^ + + xy + - 2,x - 2y. 

We want to show that P{x, y) < on (0, 1) x (0, 1). We have 

Q2{y) = + 
Oi(y) = -3 + y + 2y2, 

Qo(y) = -22/ + y'. 

First, we have Q2{y) < 1 + + (y^ — y"^) (by rising the power of —y). Now, we 
redefine Qi(y) = (— 3+l)+i/+2y^ and wc get the bound Qi(y) < 0+{y—2y)+2y^ < 
+ + (2y2 - 2/2). Hence, Qo{y) = ~ 2y + {y"^ + y^) and as the last step 
Qo{y) < + + (2t/2 - 2t/2) = 0. This shows that the inequaUty is true. 

7. Mathematica package for variational bounds 

The package Triglnt.m contains functions helpful in finding variational upper 
bounds. All functions have short explanations available (use ?Fun in Mathematica). 
This package requires Mathematica 6.0. Although all functions are written for 
triangular domains, one could triangulate any polygon and still use the pac;kagc. 

The function Trigint is equivalent to Integrate, but it is much faster for lin- 
ear combinations of trigonometric functions. It is necessary due to a very slow 
integration of some trigonometric functions in Mathematica 5.1 and above. 

(* :: Package : : *) 
BeginPackage [ " Trigint ' " ] 
Trigint : : usage = 

"Trigint is faster than Integrate for large trigonometric functions 



Equilateral : : usage= 

" Eigenfunctions of Neumann and Dirichlet Laplacian on the 
equilateral triangle with vertices (0,0), (1,0), (1/2, 
Sqrt[3]/2) : 

Equilateral I Dirichlet , Symmetric] [m,n] — l<=rrK=n 
Equilateral [ Dirichlet , Antisymmetric] ]m,n] — l<=rrKn 
Equilateral ]Neumann, Symmetric] [m,n] — 0<=iTK=n 
Equilateral [Neumann, Antisymmetric] [m,n] — 0<i=nKn 
Equilateral [ Eigenvalue] [m,n] 
Eigenfunctions of the right triangle with vertices (0,0), (1,0), (0 , Sqrt 
IS]) : 

Equilateral [ Dirichlet , Half ] [m,nj — l<=rrKn 
Equilateral [Neumann, Half J [m^n] — 0<=rrK=n 
Equilateral [ Eigenvalue , Half ] [m,n] 



Vertices : 



Square : : usage= 



Equilateral [ ] 
Equilateral [ Half] ' 



" Eig en functions of Neumann and Dirichlet L aplacian on the 
square with vertices (0,0), (1,0), (1,1), (0,1): 
Square [ Dirichlet ] [m,nj — nC>=l,n>=l 
Square [Neumann j [m,n] — m>=0,n>=0 
Square [Eigenvalue ] [m, nj 
Eigenfunctions of the right isosceles triangle with vertices (0,0), (1,0), 
(0,1): 
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29 Square [ D irichl et , H alf ] [m,n] — K^rrKn 

30 Square [Neumann f Half ] [m,n] — CK=tr<=n 

31 Vertices : 

32 Square [ Half ]" ; 
33 

34 Rayleigh : : usage="-Ra3/^e^5/^ quotient for functions on triangular domains. 



35 Many functions can be specified , each having more than one domain, hut all 
domains should be disjoint . 

36 

37 Rayleigh [fl , Til, T12, . . . , Tln,f2 , T21 , T22 , . .. , T2m, ...]"; 
38 

39 Transplant :: usage="I/mear change of coordinates . 
40 

41 Transplant [ Function , Tar get Triangle , InitialTriangle /"; 
42 

43 Limits :: usage=" Generates limits for integration over a triangle with one 
side contained in the x—axis."\ 

44 

45 Del :: usage=" Gradient with respect to x and y."; 

46 Grad :: usage=" Gradient with respect to x and y."\ 
47 

48 Area : : usage="Area of a triangle with given vertices."; 

49 Perimeter :: usage="Perimeter of a polygon with given vertices."; 

50 T : : usage="T[ a , b ] — triangle with vertices (0,0), (1,0) and (a,b)"; 
51 

52 Begin [" 'Priuate "'] 
53 

54 (* Trigint *) 

55 Trigint [ f_ , x__] := Expand [ Tint [TrigReduce [ f ] , x ] ] ; 

56 Tint [ f- , x__ , y_] i=TInt [Tint [ f ,y] ,x ] ; 



57 Trigint :: nomatch="Aro match for '1'. Integrate used."; 

58 ('* integrals *J 

59 psin [ n_] := psin [n] = Evaluate [ Integrate [y"n Sin[#l+#2 y ] , { y , #3 , # 4 } ] ] & ; 

60 pcos [ n_] : = pcos [n] = EvaIuate [ Integrate [y"n Cos[#l+#2 y ] , { y , #3 , #4}] ] & ; 

61 PP [ n_] : = pp [n] = Evaluate [ Integrate [y~n,{y,#l,#2}]]&; 

62 (* substitutions * ) 



63 sub [ x_ , a- , b-] : = { 

64 X_. x-n_. Sin[A_. + B_. x ] / ; FreeQ [ { X, A, B , n } , x] : > (X ps i n [ n ] [ A , B , a , b ] ) , 

65 X_. x"n_. Cos[A_. + B_. x j / ; FreeQ [ {X, A, B, n } ,x] : > (X pcos [ n ] [ A, B , a , b ] ) , 

66 X_. Sin[A_. + B_. x] / ;FreeQ[{X,A,B} ,x]->(X psin [0] [A,B, a ,b] ) , 

67 X_. Cos[A_. + B_. x] / ;FreeQ[{X,A,B} ,x]->(X pcos [0] [A,B, a,b] ) , 

68 X_. x-n_./;FreeQ[{X,n} ,x]: >(X pp[n] [a,b]) , 

69 X_/;FreeQ[{X} ,x]->X (b-a), 

70 X_: > (Message [Trigint : : nomatch ,X] ; Integrate [X, { x , a , b } ] ) 

n }; 

72 (* single integral *J 

73 Tint [ f_ ,{ x_ , a_ , b_}] : = ( 

74 ff=f /. Sin [X_]: >Sin [ Collect [X, { x} , Simplify ]] / . Cos [X_] : >Cos [ Collect [X,{x} , 

Simplify ] ] ; 

75 ff=Expand[ ff ] ; 

76 ff=ff+fO+fl ; 

77 ff=Replace[ff ,sub[x,a,b] ,1]; 

78 ff /.fO->0/.fl->0 

79 ) ; 
80 



81 (* eig enfunctions *) 

82 h = l; 

83 r=h/(2Sqrt [3] ) ; 

84 u=r — Global ' y ; 

85 v=Sqrt [3] /2( Global ' x-h/2) +(Global 'y-r) / 2; 

86 w=Sqrt [3]/2(h/2-Global 'x)+(Global 'y-r) /2; 

87 EqFun[f_ ,g_]:=f [Pi (-#l-#2) ( u+2r ) / (3 r ) ] g [ Pi (#l-#2) ( v-w) / (9 r ) ] + 

88 f[Pi #1 (u+2r)/(3r)]g[Pi (2#2+#l) ( v-w) / (9 r ) ] + 

89 f[Pi #2 (u+2r)/(3r)]g[Pi (-2#l-#2) ( v-w) / ( 9 r ) ] ; 
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90 Equilateral [Global ' Neumann , Global ' Symmetric] — Evaluate [ Simplify [ EqFun [ Cos , 

Cos]]]&; 

91 Equilateral [Global' Neumann , Global ' Antisymmetric] = Evaluate [ Simplify [ EqFun [ 

Cos, Sin]]] &; 

92 Equilateral [Global 'Dirichlet , Global' Symmetric] = Evaluate [Simplify [EqFun [ Sin , 

Cos]]]&; 

93 Equilateral [Global 'Dirichlet , Global' Ant isy mmet r ic] = Evaluate [ S implify [ EqFun [ 

Sin, Sin]]] &; 

94 Equilateral [Global ' Eigcn valuo] = Evaluate [4/ 2 7 ( Pi/ r ) "2(#1"2+#1 #2+#2"2)]&; 

95 Equilateral [Global ' Neumann , Global ' Half]— Equilateral [Global ' Neumann , Global ' 

Symmetric] /. Global ' x->(l-Global 'x) /2/. Global 'y->Global ' y /2// Simplify ; 

96 Equilateral [ Global ' Dirichlet , Global ' Half]= Equilateral [ Global ' Dirichlet , 

Global ' Antisymmetric ] / . Global ' x->(l- Global ' x) / 2/ . Global ' y->Global ' y/2// 
Simplify ; 

97 Equilateral [Global 'Eigenvalue , Global ' Half ] = Evaluate [ 1 / 2 7 ( Pi/ r ) "2(#1"2+#1 

#2+#2-2)]&; 

98 Equilateral []=T [1/2 , Sqrt [3] /2] ; 

99 Equilateral [ Global ' Half]=T[0 , Sqrt [ 3] ] ; 
100 

101 Square [Global ' Dirichlet] = Evaluate [Sin[#l Pi Global ' x ] Sin [#2 Pi Global'y]]&; 

102 Square [Global 'Neumann] = Evaluate[Cos[#l Pi Global ' x ] Cos[#2 Pi Global 'y]]&; 

103 Square [Global ' Eigenvalue] = Evaluate [ Pi"2(#l "2 + #2"2) ] & ; 

104 Square [ Global ' Dirichlet , Global ' Half] = Evaluate [ Sin[#l Pi Global ' x ] Sin [#2 Pi 

Global 'y]-(-l)~(#l+#2)Sin[#2 Pi Global ' x] Sin[#l Pi Global 'y]]&; 
106 Square [ Global 'Neumann, Global ' Half] = Evaluate [Cos[#l Pi Global ' x ] Cos[#2 Pi 

Global 'y] + (-l)-(#l + #2)Cos[#2 Pi Global ' x ] Cos[#l Pi Global 'y]]&; 
106 Square [Global ' Half]=T[0 , 1] ; 
107 

108 (* Transplantation 

109 LT[{ pl_ , p2_ , p3_} ,{ ql_ , q2_ , q3 _}]:= Module [ { f f } , 

110 f f [ x_] : = {x . { aa , bb} + cc , x . { dd , cc}+ f f } ; 

111 AffineTransform[{{{aa ,bb} ,{dd,oc}},{cc,ff }}]/. Solve [{ f f [pl]==ql , ff [p2]==q2 , 

ff [p3]==q3} ,{aa , bb , cc ,dd , oc , ff } ] [ [ 1 ] ] 

112 ] ; 

113 ST[p_ , q_]:=Thread[{ Global ' x , Global 'y}->LT[p ,q] [{ Global 'x , Global 'y } ] ] ; 

114 Transplant [ f _ , Tl_ , T2_] : = f / . ST [Tl , T2 ] ; 
116 

116 (* Rayleigh quotient *J 

117 Rayleigh[p__]: — Num [ p ] / Denom [ p ] ; 

118 Num[ f_ , T__List ]: = Total [ GInt [ f , #] &/@{T } ] ; 

119 Numj f_ , Longest [T__List],p__] :=Num[ f ,T]+Num[p ] ; 

120 Denom[ f_ , T __L i s t ]: = Total [ NInt [ f ,#] &/@{T } ] ; 

121 Denom [f_ ,Longest[T__List] ,P--] := Denom [ f ,T]+Denom [p ] ; 

122 GInt [ f_ ,T_]:= Trigint [Del[ f ] . Del [ f ] , Limits [T] ] ; 

123 NInt [ f_ , T_] : = Trigint [ f "2 , Limits [T] ] ; 
124 

126 Limits[{{c_,0},{d_,0},{a_, b_ }}]: = Sequence [{ Global 'y ,Min[0 ,b] ,lVIax[0 ,b] } ,{ 

Global 'x,Min[c , d] + (a-Min [ c , d ] ) Global ' y /b ,lVIax[ c , d] + (a-lVIax[ c , d ] ) Global 'y/ 

b}]; 

126 Limits [{{ c_ ,0} , { d_ , } , { a_ , b _ } , cond _}]: = Sequence® ©Refine [{ { Global ' y ,Min[0 ,b 

] ,Max[0 ,b] } ,{ Global ' x ,Min [ c , d] + ( a-Min [ c , d ] ) Global 'y/b,Max[c ,d] + (a-Max[c 
,d]) Global 'y/b}},cond] ; 

127 (* 

128 Limits [{{ c- ,0} ,{ d- ,0} ,{ a- , b -}}]: = Sequence [{ Global ' y , , b} ,{ Global ' X , c + (a—c) 

Global 'y/b , d+(a-d) Global 'y/b}]; 

129 * ) 

130 Del[ f_]: = {D[f , Global 'x] ,D[f , Global 'y]}; 

131 Grad[ f_]: = Del[f ] ; 

132 T[a_,b_]: = {{0,0},{l,0},{a,b}}; 

133 T[a_ , b_ , c_]: = {{0 ,0} ,{1 ,0} ,{a,b} ,c}; 

134 Area [{ pl_ , p2_ , p3_ }]: = Abs [ Cross [Append [ p2-pl ,0] , Append [ p3-pl ,0]]][[3]]/2// 

Simplify ; 

136 Perimeter [ 1 _]: = Total [ Sqrt [#.#]&/©( 1-RotateLeft [ 1 ] ) ]// Simplify ; 

136 
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137 End [ ] 

138 

139 EndPackage [ ] 

8. Script for triangles 

In this section we give a script that performs all the calculations from Sections [5] 
and |6] It is important to note that all operations are done symbolically, thus there 
are no numerical errors in any calculation. This allows us to use the script as a 
part of our proofs. 

The script has many comments to help the reader follow the code. The output 
of the script also contains many values representing different stages of calculations. 
The meaning of each value is explained in the comments. 

The script handles both bounds from Theorem It can be executed either 
inside Mathematica GUI or using a command line 

MathKcrncI —run " bound— _:« s cri p t .m" 

where the value of the variable "bound" should be either 1 or 2. The package 
Trigint should be in the same folder as the script. The script generates all cases 
considered in the proof. Each case corresponds to one of the test functions. For 
the second bound, the second test function is not needed, hence we used it to split 
the rectangle for the third test function. This is the only case when the algorithm 
fails to give a proof for a whole rectangle. 

Note that long lines are split into multiple lines and continuation is indented. 

1 ('* bound: 1 or 2 ^) 

2 If[! ValueQ [ bound ] , bound = 1]; 
3 

4 AppendTo [ $Path , NotebookDirectory [ ] ] ; 

5 « Trigint ' 

6 f** test functions ^ ^ ) 

7 ('* the equilateral triangle *J 

8 f[l] — EquilatcraI[Dirichlct, Symmctric][l, 1]; 

9 f[S21] = Equilateral [ Dirichlot , Symmetric ][ 1 , 2]; 

10 f[A21] = Equilateral [ Dirichlot , Antisymmetric ][ 1 , 2]; 

11 the right isosceles triangle *J 

12 g[l] = Square [ Dirichlot , Half][2, 1]; 

13 g[2] = Square [ Dirichlot , Half][3, 1]; 

14 the half of the equilateral triangle ^) 

15 h[l] = Equilateral [ Dirichlot , Half][l, 2]; 

16 h[2] = Equilateral [ Dirichlot , Half][l, 3]; 

17 (^^ find maximum of a rational function 

18 ff = (c [1] + c [2] Alpha + c[3] Alpha ~ 2 ) / ( d [ 1 ] + d[3] Alpha" 2); 

19 ff = ff /. Solve[D[ff , Alpha] = 0, Alpha] // FullSimplify; 

20 ff = ff [ [2] ] ; 

21 (^^ algorithm for solving inequalities 

22 * implementation of Section |6j * * ^ 

23 CumFun[f_, 1_] := Rest [ FoldList [ f , 0, I]]; 

24 PolyNog[P_, {x_, y_}, {dx_,dy_}] := ( ( Fold [CumFun [Min[# 1 , 0]/dy + #2 &, IVD^ 

[Max[#l, 0] &, #1] dx + #2] &, 0, Reverse [ CoefficientList [P , {x , y}]]] 
// Max) <= 0) ; 

25 

26 Do[ 

27 Print[" case; ", case]; 

28 ^ case 

29 fun[l] = Transplant [Alpha f[l] + f[S21], T[a, b], E q u i 1 a t e r a 1 [ ] ] ; 

30 (^rectangles for each bounds) 

31 r[l][l] = {{103/100, 36/100, 1, 1.45}, {0, 3/100, 0, 0.15}}; 

32 r[2][l] = {{100/100, 37/100, 1, 1.5}, {0, 9/100, 0, 0.5}}; 
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33 (** case 2* * ) 

34 fun[2] = Transplant [Alpha f[l] + f[A21], T[a, b], Equilateral []] ; 
36 r[ll[2l = {{1, 3/100, 1, 1.5}, {0, 2/10, 0, 0.3}}; 

36 (**case S**) 

37 fun [3] = Transplant [Alpha h[l] + h [ 2 ] , T[a, b], Equ 11 a t e r al [ Half ] ] ; 

38 r[l][3] = {{139/100, Infinity, 1, 2}, {0, 1, 0, 1}}; 

39 (*two different rectangles for the second bound*) 

40 r[2][3] = {{137/100, 63/100, 1, 2.5}, {0/100, 22/100, 0, 0.5}}; 

41 If [case = 2 && bound = 2, case = 3; 

42 r[2][3] = {{137/100, 63/100, 1, 2.5}, {22/100, 20/100, 0, 0.5}}]; 

43 (** case 4** ) 

44 fun [4] = Transplant [Alpha g [1] +g[2], T[a, b], Square [ Half ]] ; 
46 r[ll[4] = {{1, 39/100, 1, 2}, {1/5, 4/5, 0, 1}}; 

46 r[2][4] = {{1, 37/100, 1, 1.5}, {20/100, 22/100, 0, 0.5}}; 

47 (**case 5**) 

48 fun[5] = fun[l] + 1/2 Tr an splant [ g [ 2 ] , T[a, b], Square [ Half ]] ; 

49 r[ll[5] = {{103/100, 36/100, 1, 1.5}, {3/100, 17/100, 0, 0.25}}; 

60 r[2][5] = {{100/100, 37/100, 1, 1.5}, {9/100, 11/100, 0, 0.5}}; 

61 (**test if on the sides of the triangle**) 

62 (* output 1 — 3: three zeros if boundary conditions are met*) 

63 fun[case] /. y -> // FullSimplify // Print; 

64 fun [case] /. y -> b x/a // FullSimplify // Print; 

66 fun [case] /. y -> b (1 - t) / . x -> a - ( a - 1 ) t // FullSimplify // 
Print ; 

66 (**bounds for the eigenvalues**) 

67 frac = Simplify[b"2 Rayleigh [ fun [ case ] , T[a, b,b>0]], {b > 0, a <= 1/2}] 

// Together; 

68 (•-^•-^test if integral has a correct form, rational function**) 

69 outpui 4-riumber of the coefficients of the numerator*) 

60 (cN — CoefficientList [Numerator [ frac ] , Alpha]) // Dimensions // Print; 

61 (^output 5 : CO efficients of the denominator*) 

62 (cD = C oe fficient L is t [ Denominator [ fr ac ] , Alpha]) // Print; 

63 frac = ff /. c[i_] :> cN [ [ i ] ] /. d[i_] :> cD [ [ i ] ] ; 

64 lainbda2 = frac /b " 2 ; ('* i/pper bound*) 

66 diam = If[case !— 2, N, 1] ;(* diameter* ) 

66 R=b/(l+M+N);('*inradi«s*J 

67 lambdal[l] = Pi"2 (4/(diam"2) + ( diam " 2 ) / (M" 2 - a"2)); (*lower bound: 

Freitas * ) 

68 lambdal[2] = 8 Pi " 2/ Sqrt [ 3 ] / b ; (*lower bound: Polya*) 

69 rM = r [bound ] [ case ] [ [ 1 ] ] ; 

70 rU = r [bound ] [ case ] [ [ 2 ] ] ; 

71 bd[l] := FullSimplify [27/16/Pi-2 (lambda2 - lambdal [ 1 ] ) R-2 /. b -> Sqrt [ 

M-2 - a "2] /. a -> (M-2 - N-2 + 1) /2] ; 

72 bd[2] := FullSimplify [3 lambda2/7/ lambdal [ 2 ] /. b -> Sqrt [M-2 - a-2] /. a 

-> (M-2 - N~2 + 1) /2] ; 

73 bd2 = Numerator [ bd [ bound ]]" bound — Denominator [ bd [ bound ]]" bound ; 

74 bd3 = If [case != 2, bd2 , bd2 /. {M -> 1/N, N->M/N}] /. N -> U + M; 

76 RegionPlot [ { bd3 < / . N -> U + M, rM [ [ 1 ] ] <= M <= rM [ [ 1 ] ] + rM [ [ 2 ] ] && 
rU[[l]] <=U<=rU[[l]] +rU[[2]]}, {M, rM [ [ 3 ] ] , rM [ [ 4 ] ] } , {U, rU 
[[3]] , rU[[4]]}] // Print; 

76 (** change inequality into polynomials**) 

77 bd2 = If [case != 2, bd2 /. N -> U + M, bd2 /.M->N-M+1 /. N->1- 

u]; 

78 bd2 = bd2 /. M->M+rM[[l]] /. U->U+rU[[l]] // ExpandAll ; 

79 p= Position [ bd2 , Sqrt [ _ ? ( ! NumberQ[#] &)]];('* root expression*) 

80 PI = ReplacePart [ bd2 , 0, p] // ExpandAll ;('* poij/nomioi part*) 

81 P2 = (bd2 — PI) "2 // ExpandAll ; expression under the root*) 

82 (* output 7—8:check if we have a polynomial*) 

83 PolynomialQ [PI , {M, U}] // Print; 

84 PolynomialQ [P2, {M, U}] // Print; 
86 (**proofs**) 

86 (* output 9— 10: symbolic solution*) 

87 PolyNeg[Pl, {U, M} , {rU [ [ 2 ] ] , rM [ [ 2 ] ] } ] // Print; 

88 PolyNeg[P2 - Pl-2, {U, M} , {rU[[2]] , rM [ [ 2 ] ] } ] // Print; 
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89 , { case , 5 } ] ; 
90 

91 Exit[] 
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